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A CFD validation workshop for synthetic jets and turbulent separation control (CFD- 
VAL2004) was held in Williamsburg. Virginia in March 2004. Three cases were inves- 
tigated: synthetic jet into quiescent air. synthetic jet into a turbulent boundary layer 
crossflow, and flow over a hump model with no-flow-control, steady suction, and oscilla- 
tory control. This paper is a summary of the CFD results from the workshop. Although 
some detailed results are shown, mostly a broad viewpoint is taken, and the CFD state- 
of-the-art for predicting these types of flows is evaluated from a general point of view. 
Overall, for synthetic jets, CFD can only qualitatively predict the flow physics, but there 
is some uncertainty regarding how to best model the unsteady boundary conditions from 
the experiment consistently. As a result, there is wide variation among CFD results. For 
the hump flow. CFD as a whole is capable of predicting many of the particulars of this 
flow provided that tunnel blockage is accounted for, but the length of the separated region 
compared to experimental results is consistently overpredicted. 

I. Introduction 

The study of separation control through the use of blowing, suction, and synt hetic jets has been an act ive 
field of research for some time now. 1 * Many experimental and theoretical papers (only a few of which are 
referenced here) have been published on the subject. 2-13 The review article by Greenblatt and Wvgnanski 7 
lists many additional references. There has also been a recent increase of attention paid to the field by the 
American Institute of Aeronautics and Astronautics (AIAA), which held its first Flow Control Conference 
in June 2002 in St. Louis, Missouri, and now sponsors it as a continuing biennial series. 

Computational methods have also been used extensively for computing synthetic jet flows. Only a few 
examples are referenced here. 8 13 However, to our knowledge most CFD efforts have been somewhat isolated, 
so it has been difficult to make an assessment of the state-of-the-art as a whole. In order to more broadly 
assess the current capabilities of different classes of computational methodologies, the time appeared to be 
ripe to hold a workshop for which many participants would compute the same flow control test cases. To 
this end, the CFD Validation of Synthetic Jets and Turbulent Separation Control (CFDVAL2004) workshop 
was held in Williamsburg, Virginia in March 2004. 14 Three different test cases, all of which were carried out 
experimentally at NASA Langley Research Center, exercised various aspects related to the flow r physics of 
separation control. 
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This paper is a summary report of the overall CFD results from the CFDVAL2004 workshop. It does 
not, direct ly address technical issues associated with the experiments themselves, except as they relate to the 
CFD modeling of the flows. Technical papers describing the experiments are also available. 15-17 This paper 
is organized as follows. First, a brief description of the workshop is given, including an overall description 
of the three test cases. Then, a summary of results from each of the three test cases is presented. Finally, 
the paper closes with summary and conclusions. 

II. The CFDVAL2004 Workshop 

The CFDVAL2004 workshop was structured similarly to the series of Refined Turbulence Modeling work- 
shops sponsored by the European Research Community on Flow, Turbulence, and Combustion (ER.COF- 
TAC). The 10th such workshop was held in October 2002 in Poitiers, France. Although not formally a 
part of the European series, CFDVAL2004 was held in association with the ERCOFTAO organization, and 
representative members from that organization contributed to the workshop. CFDYAL2004 was also asso- 
ciated with the Air Force Office of Scientific Research (AFOSR), the International Association of Hydraulic 
Engineering and Research (IAHR), QNET-CFD, and the National Institute of Aerospace (NIA). 

There were 75 attendees at the workshop. Seven countries were represented, including the United States, 
France, Italy, Germany, Japan. United Kingdom, and Switzerland. Most of the workshop participants came 
from universities, but several companies and public sector research laboratories were represented as well. 

The purpose of CFDYAL2004 was to bring together an international group of computational fluid dynam- 
ics practitioners to assess the current capabilities of different classes of turbulent flow solution methodologies 
to predict flow fields induced by synthetic jets and separation control geometries. The experimental data 
for the workshop were taken specifically with CFD validation in mind. Therefore, a great deal of effort, was 
expended to attempt to document most of the relevant boundary condit ions. The 3-day workshop gave time 
for each participant to discuss their method and results. Each test case was concluded by a summary talk 
by the organizers, for which all of the contributed results were compiled. A panel discussion, whose focus 
was quality and trust for synthetic jet CFD and experiment , was also included in the workshop. 

To encourage broad participation and to determine the general state-of-the-art, the decision was made 
early not to dictate particular boundary condit ions, grids, or method of solution. For example, although 
the experimental velocity was given as a function of time near the centers of t he jet exits for cases 1 and 2, 
the participants were given the choice as to the precise boundary condition used to try to match it. While 
this strategy to broaden workshop participation was successful - solution methods ranged from reduced- 
order models through Reynolds-averaged Na.vier-Stokes (RANS), to large-eddv simulation (LES), and direct 
numerical simulation (DNS) - it also had an unavoidable downside. Differences in grids and boundary 
conditions represented an additional source of uncertainty when attempting to compare CFD results with 
each other. 

The three test cases were chosen to represent different aspects of flow control physics. In test case 1 
(synthetic jet. into quiescent air), flow passed in and out of a slot (1.27 mm wide by 35.5(5 mm long), which 
was located on the floor of an enclosed box 0.(51 m per side. The jet was driven bv a side-mounted circular 
piezo-electric diaphragm inside the cavity chamber beneath the floor. The frequency was approximately 445 
Hz, and the maximum velocity out of the slot was approximately 25-30 m/s. T his case was considered to be 
nominally two-dimensional at the center plane of the slot. For case 1, the participants were asked to supply 
long-time-average jet width and velocity profiles at several locations. Also, phase- aver aged velocities and 
turbulence quantities were requested at 8 different phases at several locations. Finally, time-history values 
of velocity at three points in the flow, as well as specific line contour plots of phase- averaged velocities were 
also requested. 

In test case 2 (synthetic jet in a crossflow), flow passed in and out of a circular orifice (5.35 mm in 
diameter. The orifice was located on the floor of a wind tunnel splitter plate with a turbulent boundary 
layer at M=0.1 and approximate boundary layer thickness of 21 mm. The jet was driven elect ro-mechanically 
by a bot tom-mounted square-shaped rigid piston mounted on an elastic membrane inside the cavity chamber 
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beneath the splitter plate. The cavity was approximately 1.7 mm deep and the piston moved approximately 
±0.77 mm. The frequency was 150 Hz, and the maximum velocity out of the orifice was approximately 
3 .3(7^. For case 2. both long-time- average jet width and velocities as well as phase-averaged velocities and 
turbulence quantities (9 phases) were requested at several locations. Furthermore, participants were asked 
to supply time-history values of velocity at three points in the flow, as well as specific line contour plots of 
bot h long- t ime- average and phase- averaged u- velocity. 

In test case 3 (flow over a hump model), turbulent flow at M— 0.1 passed over a hump of chord 420 mm 
mounted on the floor of a wind tunnel splitter plate. The hump used end-plates at both sides (the model 
had a span of 584.2 mm between the plates), and the flow was nominally two-dimensional at its center 
plane. The hump had a slot near 65% chord, near where separation naturally occurred. This case had two 
mandatory conditions: no-flow-control (no forced flow through the slot), and steady suction (m = 0.01518 
kg/s). There was also one optional condition (the experiment was not completed in time for the workshop) 
of oscillatory (synthetic) jet control. In this case, the jet was driven electro- mechanically by a bottom- 
mounted rectangular-shaped rigid piston mounted on an elastic membrane deep inside the cavity chamber. 
The frequency was 138.5 Hz, and the peak velocity out of the slot was approximately 27 m/s. For the case 3 
no- flow-control and steady-suction conditions, participants were asked to supply Cp and Cf along the hump 
surface, velocity and turbulence profiles at several locations, and streamline contour plots. For the optional 
synt hetic jet condition, only long-time- average Cp along the hump surface was requested. 

There was a great deal of CFD data submitted to this workshop. It is clearly not possible to show 
most of the results here. Summary plots as well as several specific detailed plots are shown that we believe 
to be representative of the solutions as a whole, or which serve to illustrate specific points or differences. 
Additional detailed plots of the CFD results are available from the CFDVAL2004 website. 14 It should also 
be noted that workshop participants were allowed to correct, resubmit, or withdraw their submissions during 
the month following the workshop; both this paper summary and the plots on the website contain updated 
information as of May 2004. 

III. Case 1: Synthetic Jet into Quiescent Air 

Tabulations of the submissions for Case 1 are given in Tables 1 and 2. There were 8 contributors who ran 
25 separate cases. There was one LES submission, one reduced-order model submission, one 2-D blended 
RANS-LES submission, and several laminar Navier-Stokes submissions (the 3-D laminar runs perhaps can 
also be characterized as under- resolved DNS); the others were unsteady RANS (UR A NS). Most of the runs 
were computed in 2-D, but there were also a few 3-D computations. Note that none of the 3-D comput ations 
modeled the actual shape of the the cavity, including the circular diaphragm; instead, they were all computed 
using periodicity in the direction aligned with the slot 's long axis. Six of the contributors modeled the cavity 
(or some approximation of the cavity) and two did not model any cavity. Of those who modeled the cavity, 
ONER A-flu3m and UKY-ghost applied a time- varying velocity on the side of the cavity where the diaphragm 
was located, derived from diaphragm-center displacement data. WASHU-wind and NASA-t lns3d also applied 
a similar boundary condition there, but based it on best matching of the data at the slot exit. GWU-vicar3d 
used a different cavity shape with time- varying velocity boundary condition applied at the bottom wall. 
NCAT-quasld modeled the motion of the actuator with a quasi- 1-D model. 

All of the results were calculated with structured grids. Only two contributors (UKY-ghost and NASA- 
tlns3d) ran with the same grid, U 2-D Structured Grid #1". This grid, along with other structured and 
unstructured grids, were made available on the website several months prior to the workshop. Descriptions 
of these grids can be found at the website. 14 

In this section, w 7 e compare results for case 1 with particle image velocimetrv (PIV) and hot wire data. 
It should be noted that other experimental data were later taken using a different piezo-electric diaphragm 
(these devices are prone to occasional failure). This later data also include laser doppler velocimetrv (LDY) 
measurements in addition to PIV and hot wire, but they w r ere also at slightly different conditions (higher 
maximum outflow 7 velocity) than the data used for the workshop. The later data are not shown here, but 
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are reported in Yao et al. 15 A sketch showing the locations where most, of the results from the computations 
were supplied is given in Fig. 1, and the definition of jet width used for this case is given in Fig. 2. 

Time histories of v- velocities near t he center of the slot, exit are shown in Fig. ■>. There is clearly a 
significant variation between the two measurement techniques (PIY and hotwire), so it, is difficult to make 
a confident assessment of the CFD results. Both experimental measurements indicated approximately a 
200° phase difference between peak and trough, while many of the CFD methods yielded 180°. Some of 
the differences in particular phase- averaged results between CFD and experiment could be a result of not 
matching the same phase. 

Fig. 4 shows long-time-average v- velocity profiles along the centerline, x=0. Again, a significant differ- 
ence can be seen between the two measurement techniques, especially very near the slot for v less than 
approximately 2 mni. There were also significant variations between the CFD results. Most of them showed 
similar t rends to the experiment, although NCAT-quasld and WARWICK-neat indicated a more rapid drop 
in velocity for y > IU mm compared to the other results. 

Long-time-average results along two y=constant. lines (y=().l mm and y=4 nun, respectively) are shown in 
Figs. 5 and (). Just above the slot exit, at y=0.1 mm, all results except POlT-saturne greatly underpredicted 
the average jet magnitude. However, beginning near 1-2 mm and above, all the CFD results agreed very well 
with the PIV experiment. An example of this better agreement is shown at y—4 mm in Fig. (). There was 
relatively little spread among the CFD results at this location. Although not shown, the variation between 
the CFD results increased at higher y-stations. 

The long-time-average jet width is shown in Fig. 7. Most CFD results tended to somewhat overpredict, 
the width. The OXER A-flu3m laminar result was similar to the other methods near to the wall, but then 
predicted too much spreading of the jet past y=4 mm. It, was more than twice as wide as the experimental 
results near v=8 nun. Although not shown, the CFD results had significantly greater variation for v > 8 
mm. 

For the sake of brevity, only one representative phase- averaged result is shown here. Fig. 8 shows v- 
velocitv profiles at v=4 nun, at phase=135°. Compared with the long-time-average results shown earlier, 
this typical phase- averaged result illustrates that there was a great, deal more CFD variation at, this level of 
detail. Also, the comparison was not as good with experimental data. 

Several participants investigated the effect, of time step and grid size on the solution. Although not 
shown, time step had very little effect, for any of the solutions. The effect, of grid was relatively small near 
the wall, but it could be larger away from the wall. This is illustrated in Fig. 9, which shows long-time- 
average v- velocity along the centerline. UKY-ghosfc and POIT-saturne showed only minor variation, but the 
3-D GWU- vicarSd results were affected by grid size and its spanwise extent. NASA-tlns3d fine and medium 
grids gave nearly identical results, but the coarse grid was different,. 

The turbulence model s effect on the CFD results was found to be fairly significant for this case. Results 
are shown in Fig. 10. For POIT-saturne, the k-s model compared better with experimental data than the 
Reynolds stress model (RSM) for this quantity. For WARWICK-neat, three different turbulence models 
were essentially the same for y < 4 mm, but then they behaved very differently beyond that. None of these 
agreed particularly well with experimental data. For WASHU-wind, Spalart-Allmaras (SA) and Mentor's 
shear stress transport (SST) were fairly similar, especially for y < 7 mm, but the 2-D blended SST-LES 
method generally overpredicted the average velocity everywhere. For NASA-tlns3d, SA and SST showed 
similar variation to that seen by WASHU-wind: the two models were close near the wall, then SST predicted 
somewhat larger velocities than SA further from the wall. 

Although not shown, the computed turbulence quantities (e.g., u'v') did not compare well at all with 
experimental data, and there was very large variation exhibited among CFD results. Also, there is some 
question whether the measured data is entirely a result of turbulence, or a result of unsteady “flapping" 
of the time-dependent jet. It, was generally felt, by the participants of the workshop that this case 1 flow 7 
is probably laminar in the beginning, then transitional further into the flow field. Therefore the use of a 
standard R ANS turbulence model may not be entirely warranted. 

As mentioned earlier, in spite of attempting to match the v-velocity phase at t he slot outflow 7 (see Fig. 3), 
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many of the participants had difficulty because the experimental variation was not sinusoidal (18U° peak- 
topeak). As a result, at certain phases they found that the computed jet could be at a different position 
than in the experiment. This is illustrated by example flowfield contours at phase=135°, shown in Fig. 11. 
Similar to most other results (not shown), the NASA-t.lus3d result at this phase had the extent of its peak 
jet velocity lower (near 4 mm) than seen in the experimental data (near 7 mm). 

For this case, the 2-D laminar 0NERA-flu3m result predicted too much jet spreading. However, the 2-D 
laminar NCAT-quasld, 3-D laminar, and 3-D LES results were consistent with other CFD methods in this 
regard. Different turbulence models could have significant effects, especially further into the flow field: the 
SA, SST, and POIT-saturne-ke models seemed to be of similar quality, and generally agreed best with the 
data among the URANS models. Surprisingly, there did not appear to be a distinct advantage to modeling 
the cavity for this case, as opposed to specifying a jet boundary condition at the surface. Regarding grid 
effect , 2-D grid levels of about 60.000 points (including cavity) were found to be sufficient to capture the flow 
physics, whereas 16,000 points yielded noticeable differences in the solution. Among 3-D results, changing 
from 460.000 to 700,000 points (or increasing the span wise periodic extent) could significantly alter the 
solution far into the flow field. There was very little effect of changing time step, even with as low as 72 
steps per period. 


IV. Case 2: Synthetic Jet in a Crossflow 

Tabulations of the submissions for Case 2 are given in Tables 3 and 4. There were 5 contributors who ran 
10 separate cases. There was one LES submission: the others used URANS. All runs were (necessarily) 3-D. 
All methods were second order in space and time. All submissions modeled the cavity, except for CIRA-zen 
which specified a time-varying profile at the orifice exit. However, both USTO-rans and NASA-fun3d altered 
the cavity shape. In the latter case, the bottom wall was made flat while keeping the volume the same as 
the original (stepped) geometry. Two of the contributors modeled the full plane, while the others modeled 
a half plane with symmetry imposed at the center plane. Each contributor used a different grid: four used 
structured grids and one used unstructured. The “Structured Grid #1" and "Unstructured Grid #F* were 
made available on the website several months prior to the workshop. Descriptions of these grids can be 
found at the website. 14 The four contributors who modeled the cavity applied a time-dependent velocity 
specification at the cavity bottom wall: V — \ cos(27r/f). Each adjusted U to achieve what was considered 
to be a reasonable match of the experimental w- velocity variation at the orifice exit . 

In this paper, case 2 results are quantitatively compared only to LDV data. A sketch showing some of 
the LDY data measurement locations is given in Fig. 12. PIY data was also acquired at several planes. 
The reader interested in quantitatively comparing experimental results using the two different techniques is 
referred to Schaeffler and Jenkins. 

Time histories of velocities near the center of the orifice exit are shown in Fig. 13. The u- velocity 
component was captured well by three of the codes. CIRA-zen showed a higher value and NASA-fun 3d 
showed a lower value than experimental results during the expulsion part of the cycle. All codes predicted 
a zero or near-zero cross-stream v-component of velocity, whereas the experimental data indicated a very 
large v-component during expulsion. The cause of this high experimental value is not known. Finally, all 
codes captured the w- velocity reasonably well (with some variation), although no one replicated the "dip" 
near phase=160°. The cause for this dip in the experimental result is also not known. The cavity under 
the orifice was extremely shallow (nominally 1.7 mm deep neutral position with tunnel on), and the piston 
moved ±0.77 mm up and down from this position. The volume in the cavity more than doubled from 
roughly 10,000 to 26,000 cubic mm during the cycle. It is questionable whether a simple non-moving- wall 
transpiration boundary condition, used by all of the CFD met hods, is a very accurate model for this large 
a range of piston motion. 

Fig. 14 shows long-time-average u- velocity profiles on the center plane at three downstream stations. 
Overall, there was fair agreement between CFD arid experiment, but there was also a significant spread 
among the CFD results. Fig. 15 shows w- velocity profiles along one spanwise line and two stream wise lines. 
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Along t lie spanwise line, there were no experimental LDY data, and there was a significant spread in the 
CFD results. Along the streamwise lines, the 5 CFD results again exhibited significant spread; two of the 
CFD methods exhibited reasonably good agreement with experimental results along the z=10mm, y=U line 

(Fig- 15c). 

Two of the participants examined the effect of grid size on the solution. The results can be briefly 
summarized in a representative figure, Fig. lb. In general, the fine (3.9 million cells) and coarser (0.49 
million cells) grids used in NASA-cfl3d did not exhibit a dramatic influence in the solution. However, for 
NASA-fun3d the fine (0.2(5 million nodes) and coarser (4(5,000 nodes) grids did have a significant difference. 
For NASA-fun3d, the finer grid produced a negative streamwise velocity during expulsion, but the coarser 
grid did not. 

Two of the participants investigated the effects of different turbulence models on the solution. Again, 
only one representative figure is shown here for the sake of brevity, Fig. 17. In NASA-cfl3d, three different 
models resulted in minor differences from each other, and none was clearly better in comparison with the 
experiment al data. The two models used bv USTO-rans also showed only minor differences from each other, 
and they were clearly different in character from both the NASA-cfl3d results and the experimental dat a. In 
other words, the turbulence models themselves did not have as much of an impact as the combined effect of 
different codes, grids (in number of points and full plane vs. half plane), and other solution variations (such 
as alteration of t he cavity shape by USTO-rans). 

Only a few phase-averaged results are shown here. Fig. 18 shows u- velocity along a line ID downstream 
at three different phases. All CFD results were in good agreement with each other and with experimental 
results at phase=0° (at this time the orifice is still in its suction phase, so there is little to no dynamic flow 
occurring ID downstream). However, at the two later phases the influence of the expulsion part of the cycle 
is felt at this location. The CFD results captured the influence, generally following the trends seen in the 
experiment, but there was a significant spread in the solutions. 

Phase- averaged turbulence quantities from the CFD solutions generally exhibited similar trends to each 
other, but again with a significant spread. There were usually larger discrepancies between CFD and exper- 
iment than the discrepancies seen between turbulence profiles. An example is shown in Fig. 19. 

Finally, a sample comparison between experiment and CFD over the entire plane ID downstream is shown 
in Fig. 20. This figure shows u-velocity contours at phase=120°. It is intended as a qualitative comparison 
only. Also, only one CFD result is shown for brevity: other CFD results were qualitatively similar. Although 
details were different, this figure indicates that CFD predicted the overall structure of the dynamic flow 
resulting from the expulsion part of the cycle as it passes downstream. 

For this case study, LES vs. URANS did not appear to be a significant factor. In fact, overall (in 
a subjective sense) the two methods that tended to yield the most similar long-time-average and phase- 
averaged mean-flow solutions were ONER.A-flu3m (LES) and NASA-cfl3d (URANS). In general, no one 
method, algorithm, or turbulence model stood out as being the best methodology for capturing the physics 
of this flow'. Grid levels of about a half a million points or more appeared to be sufficient to capture the 
general character of the unsteady flow sufficiently well, whereas grids coarser than this could change the 
solution significantly. None of the participants submitted results examining the effect of time step; however, 
NASA-cfl3d performed computations that were not submitted, varying between 720 and 1440 steps per cycle, 
and found very little effect. 

V. Case 3: Flow over a Hump Model (Separation Control) 

Tabulations of the submissions for Case 3 are given in Tables 5 and 6. There were 13 contributors 
who ran 5(5 separate eases. Methods were mostly RANS/URANS, but there was also one DNS result and 
several blended RANS-LES results. Most of the runs were computed in 2-D, but there were also several 
3-D computations. Some of the 3-D computations were periodic in the spanwise direction, and some were 
half-plane (either flat in viscid side wall or actual side plate geometry modeled viscously). Most of the 
contributors modeled the cavity, but several applied boundary conditions directly on the hump surface. 
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Results were rom puled on both structured and unstructured grids. Many contributors used the 2-D grids 
supplied oil the website (for example, “2-D Structured Cirid #1"). Descriptions of these grids can be found 
there. 14 

Case 3 is different from cases 1 and 2 in that most of the results were not time-dependent . Only the 
optional oscillatory-control case was unsteady, and for that case only long-time-average Cp's were requested 
from the participants (8 of the contributors computed this, running 12 separate cases). However, when there 
is more extensive time-dependent data available in the future, the oscillatory-control experiment may serve 
as a good test case for time-accurate validation exercises. A sketch showing the locations where most of the 
results from the computations were supplied is given in Fig. 21. Computed pressures are compared with 
experimental data from surface- mounted pressure taps, and velocity and turbulence profiles are compared 
with PIY 7 experimental data. Experimental details are reported in Creenblatt et al. 1 ' 

A summary of all of the CFD results for the no- flow-control condition is shown in Fig. 22. (Note 
that separation and reattachment, locations are not displayed in the figures for those methods that did not 
report skin friction.) As a whole, most CFD results missed the pressure levels over the hump between 
0.2 < x/e < 0.6, and also predicted higher pressures than experimental results in the separated region 
upstream of x/c=l. As will be shown below, the missed pressures were likely due to blockage effects caused 
by the side plates in the experiment. The separation location (not known precisely for t he experiment) 
was predicted reasonably well by most CFD methods, and the reattacliinent location was predicted (for the 
most part) significantly downstream of the experimental location of x/c = 1.11. Most models of similar 
type generally behaved similarly. Two exceptions were CIRA-zen-sst-kprof-3, which predicted separation 
later and reattachnient earlier than other results with the SST model, and CIRA-zen-ke. which predicted 
reattachment further downstream than other results using a k-s model. CIRA-zen also reported unusually 
low skin friction levels for all of its results. 

A summary of all of the CFD results for the steady suction condition is shown in Fig. 23. Again, most 
CFD results overpredicted the pressure levels over the hump between 0.2 < x/c < 0.6. Also, the pressures in 
the separated region did not exhibit the same levels or shape as experimental results. The separation location 
(again not known precisely for the experiment ) was predicted to be slightly upstream of the experimental 
location, and the reattachment location was predicted to be significantly downstream of the experimental 
location of x/c = 0.94 by all CFD methods. For the suction condition, t here was generally a broader range 
exhibited in the results for a given turbulence model than for the no-flow-control condition. Notable results 
that showed the largest differences from other turbulence models of the same type were: CTR-fluent-sa-2, 
CTR-fluent-sa-1, NASA-rans-sst-weno-5, and US-fluent-ke. 

The two figures, Figs. 22 and 23, are also instructive from the point of view of CFD uncertainty. The 
surface pressure coefficient plots illustrate the range of variation of CFD for these cases. This variation is 
caused by use of different grids, codes, turbulence models, and boundary conditions (including 2-D vs. 3-D 
modeling). In particular, it should be noted that the downstream boundary condition in internal flow’ cases 
such as this can affect the flow to some degree: by setting different back pressures, one can achieve different 
results. Individual differences in this boundary condition likely account for some of the CFD variation. 

The tendency for most of the CFD results to miss the pressure levels over the front half of the hump was 
believed to be due to blockage effects in the experiment. Two pieces of evidence support this conjecture. First , 
the 3-D AZ-cobalt runs that accounted for the side plate physical shape (modeled as a viscous wall) resulted in 
improved levels there. Second, the 2-D NASA-efl3d run that used an altered top wall shape to approximately 
account for the side-plate blockage (using area ratio of plate cross-section to tunnel cross section) resulted in 
excellent agreement of pressures over the front part of the hump compared with experimental data. These 
results are shown in Fig. 24. The grid used in the NASA-cfl3d case (2-D Structured #5) is described on the 
CFDVAL2004 website; this grid was created after the workshop. It should be noted from this figure that 
the 3-D AZ-cobalt runs also yielded good results in the separated region, whereas the 2-D NASA-cfl3d runs 
did not. This suggests that there may be 3-D separated structures along the back end of the side plates that 
further constrict the flow’ (additional blockage) in the separated region. 

Although not shown, most of the contributed results w’ere in fairly good agreement with experimental data 
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(for bot h velocity profiles and t urbulent shear st ress profiles) just upstream of the slot (especially for the no- 
flow-control condition). And the velocity profiles inside the separated region were also predicted reasonably 
well, on the whole, for both no-flow-control and suction conditions. However, because the CFD results 
reattached too late, most velocity profiles compared poorly with experimental data at stations downstream 
of reattachment. The two exceptions to this were AZ-cobalt-des-l-3d and META-cfd++lns-3d for the no- 
flow-control condition. Profiles using these methods (along with five other typical representative solutions) 
at x/c = 1.2 are shown in Fig. 25. However, all results were poor downstream of reattachment, for the suction 
condition. Fourteen representative results are shown at x/c =1.0 in Fig. 20, 

A possible reason for reattachment being predicted too late is that most of the current models and methods 
predicted turbulent shear stress to be too small in magnitude in the separated region. As an illustration 
of this tendency, Fig. 27 shows bot h u- velocity and turbulent shear stress for the suction condition at the 
location x/c = 0.8, using several representative solutions. The velocity profiles exhibited a variation in 
shape, especially near the center of the profile, but the overall agreement with experimental results was 
fairly good for the BANS methods. However, the magnitude of the turbulent shear stress was seriously 
under predicted. (The turbulent shear stress data had a maximum error of as much as 20%, based on a 
pseudo-empirical uncertainty analysis. Nonetheless, the results still indicated a large difference between 
CFD and experiment.) One of the anomalies evident in the submissions was that the US-fluent results 
consistently yielded significantly larger -in-magnitude turbulent shear stress levels, in better agreement with 
experimental data, even though they used many of the same turbulence models as other codes. In any case, 
the US-fluent, solutions still predicted reattach ment too far downstream, similar to the other submissions. 

The DNS and blended RANS-LFS solutions (DFS. LNS, FSM) were generally on a par with the RANS 
methods; they, too, predicted reattachinent too far downstream. However, DNS, LES, and 3-D blended 
RANS-LFS models involve typically finer grids than RANS methods, and require very lengthy run times. It 
is possible that future improvements may result from additional efforts in these areas. As shown earlier in 
Fig. 25, a few of the blended R ANS-LES models appeared to do a bett er job predicting velocity profiles near 
and downstream of reattachment than the R ANS methods for the no-flow-control condition, but they did 
not improve results for the suction condition. Also, it is worth noting that AZ-cobalt-des predicted resolved 
turbulent shear stress levels downstream of reattachinent (for both no-flow-control and suction conditions) 
that were much higher in magnitude than most R ANS results. 

Many of the contributors investigated the effect of grid on the solution, for the no-flow-control and 
suction conditions. For 2-D, little influence was noted (all used a minimum of about 50,000 grid points). 
The influence of turbulence model has been shown earlier in Figs. 22 and 23. Generally, there were some 
differences between the models as expected, but all performed poorly when it, came to predicting reattachment 
location. No one model stood out as performing consistently better for this case. A few contributors ran 
without modeling the cavity. In general this did not appear to either help or hurt, even for the suction ca.se. 
CIRA-zen found that applying oblique suction improved results over applying normal suction. 

A summary of all of the CFD results for the (optional) oscillatory control case is shown in Fig. 28. There 
was variation among the CFD results, but in general they all exhibited similar trends. The results whose 
separated region pressure recovery occurred furthest upstream (in better agreement with experimental data) 
were UAZ-cfl3d, BOEING-overflow, and UMD-rans. NASA-fun2d investigated the effect of time step on this 
solution, and found very little effect. 

VI. Summary and Conclusions 

Time-dependent flows involving unsteady flow control are difficult not only to compute with CFD, but 
also to measure experimentally. The CFDVAL20U4 workshop in many ways ended up being as much of 
a workshop on experimental issues as on CFD issues. From the two cases involving synthetic jets (case 1 
and case 2), we learned that there are still some inconsistencies between different measurement techniques 
for certain aspects of these types of time-dependent flows. These differences were openly reported at the 
workshop, in order to foster honest and open dialog and to help encourage collaborative efforts to work 


8 of 31 


American Inst it ute of Aeronautics and Astronautics Paper 2004-2217 



toward improving future results. .Most importantly, the workshop proved to be an ideal setting for providing 
a quality technical interchange between (TO and experimental scientists, where each could learn about the 
capabilities and limitations inherent in the processes and tools of the other. 

t his paper focused on providing a summary of the CFD results from the workshop. The workshop 
itself was designed to encourage wide participation and to establish a baseline for the state-of-the-art in 
the area of synthetic jet validation and turbulent separation control. Specific grids, boundary conditions, 
and methods of solution were not dictated to the participants. As a result , CFD variation was fairly large. 
However, trends could still be noted, including overall capabilities and shortcomings. Many of the workshop 
participants suggested that, now that we have established a baseline for these cases, future CFD efforts 
should attempt to minimize known sources of CFD variation, especially regarding boundary conditions. It 
was also suggested that future experiments focus additional effort documenting the time-dependent boundary 
conditions, especially at and near the exit plane of the jet/suction slot or orifice. 

The bottom line from the CFD results ('an he summed up as follows: no one CFD technique excelled 
above others, and there was wide variation (especially for time-dependent results), and only qualitative 
agreement with experimental data. In other words, the "state-of-the-art** C FD methods of today are not 
fulh adequate to consistently and accurately predict these types of flows. What is today's "state-of-the-art"? 
Most of t he contributors used FRANS for the time-dependent cases and RANS to compute the steady cases. 
A wide variety of turbulence models were employed. Most CFD methods were second order in space and 
time. Some contributors used higher order, but there did not appear to be any obvious benefits from doing 
so for these cases. Most contributors utilized transpiration boundary conditions on a fixed grid to simulate 
moving walls. The few blended RANS-LES, LES, and DNS solutions showed merit, and were of similar 
quality to the RANS and FRANS solutions. We believe that with today's increased computer power and 
improved algorithms, these advanced methods are currently coming of age and represent the wav much of 
CFD will be performed in the future. However, as yet they showed no clear benefits over RANS/URANS in 
the sense of providing consistently better results for the workshop cases. 

Case 1 (synt hetic jet into quiescent air) was a difficult experiment to simulate. The flow field was probably 
mostly laminar or transitional, so it was unclear how best to simulate it. The piezo-electric driver used in the 
experiment was difficult to model using CFD because it moved like a drum (and possibly with multiple shape 
modes), and not like a rigid piston. The experiment indicated a deviation from periodicity in the velocity 
near the jet exit that was not completely simulated or captured by the CFD results, and caused specific 
phase results to be mis-aligned. For example, at phase— 1 35° in the experiment , the jet extended higher into 
the field than at the same phase as defined by most of the CFD runs. Long-time-average results showed 
deviations from experimental data very near to the slot as well as far from the slot. In the intermediate range 
(4-8 mm from the slot), the results were generally better. CFD yielded significantly larger variations in the 
phase- averaged quantities. Except for a few cases, average jet width was usually somewhat overpredicted 
by CFD in the near field. In the FRANS simulations, turbulence models were found to have a significant 
effect for this case. In general, the A*-f, SA, and SST models performed better compared to experimental 
data than other models tested. 


Case 2 (synthetic jet in a crossflow) was the least computed of the three cases, probably owing to its being 
both time-dependent and (necessarily) 3-D. Similar to case 1, the time-dependent experimental velocities 
measured near the orifice exit exhibited anomalies not captured or modeled by CFD. In particular, the 
experiment exhibited a dip in w-velocity on the downstroke and a very large spanwise velocity component 
during the expulsion part of the cycle that none of the CFD methods accounted for. In spite of this, 
reasonably good qualitative results were obtained compared to experimental results, but there were significant 
CFD variations in both the long-time-average and the phase- averaged results. No one method, algorithm, 
or turbulence model stood out as being the best methodology. Even LES vs. FRANS did not appear to be 
a significant factor. Generally, different turbulence models for FRANS did not have as much of an impact 
on this case as different grids, codes, and other solution variations. 

Case 3 (flow over a hump model) had the greatest number of workshop participants, since the required 
cases were nominally steady and 2-D. However, it was discovered after the workshop that the side plates used 
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in the tunnel caused blockage that, if not modeled, resulted in relatively minor (but noticeable) overpredict ion 
of the pressures over most of the hump. Also, based on some 3-D results that modeled the side plates viscously 
and obtained good pressure comparisons in the separated region, it is likely that flow structures near the back 
of the plates constrict the flow even further in that region. In other words, in spite of this flow appearing to 
be relatively simple, computing the wall pressures accurately can require full 3-D modeling or else some sort 
of blockage corrections. Overall, case 3 was not a severe test from the point of view of predicting separation. 
But GFD results were deficient in another important regard: they consistently predicted the reattaehment 
location to be significantly further downstream than the location documented in the experiment. This same 
behavior occurred regardless of turbulence model or method; even a DNS computation predicted too long 
of a separation bubble. Inside the bubble itself, most computations predicted velocity profiles in reasonably 
good agreement with experimental data, but turbulent shear stresses were underpredicted in magnitude. 
This underpredict ion is consistent with delayed reattaehment, because it indicates reduced turbulent mixing 
inside the separated region compared with experimental data. 

In conclusion, the CFDVAL2UU4 workshop has established a benchmark for three different flows involving 
synthetic jets and turbulent separation control. Although the current, state-of-the-art CFD methods are 
deficient in being able to consistently and accurately predict these flows, two areas have been identified as 
key to improvement . 

• First, for synthetic jets, it is important to employ consistent- boundary conditions when multiple CFD 
methods are compared. Only then will we be able to isolate and work to correct deficiencies in the CFD 
models and methods. Implicit in this need is the requirement that experiments be used to document 
extremely detailed and accurate time-dependent flow field variables at and near the slot/orifice exits. 

• Second, for the hump model case, turbulence models (for HANS) or other methods such as LES or 
blended RANS-LES need to be improved/developed/calibrated to increase the mixing in the separated 
region and bring about earlier reattaehment and flow recovery. 
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Table 1. Case 1 Summary of Submissions 


Label 

Organization 

Authors 

Method 

ONERA-flu3m-lrs-3d 
0NERA-flu3ni-lam 
ONER A-flififin-sa 

ONERA 

M ary 

LES, 3-D 
Laminar N-S 
URANS, SA model 

UK Y-ghost-sst 
UKY-gliost-sst(fine) 

U. Kentueky 

Huang, Chen, Katam, 
Pari ini, he Beau, Huang 

URANS, SST model 
URANS, SST model, hue grid 

( i \\ l -virariid-IU^liiK') 
UWU-virarMd-.'ld 
C;VVU-vicai:kl-:id('1.. r Kl) 

GYVU 

Rupesh, Ravi, Mittal, 
R.aju, (Dallas, Cattafesta 

Laminar N-S, fine grid, 3-D 
Laminar N-S, 3-D 
Laminar N-S, larger 
domain, 3-D 

N( 'AT-quas 1 cl + rails 

NC A&T State V. k 
NASA LaRC 

Vania lee v , ( ■ ar pe n ter , 
Vatsa 

Reduced-order model 
in slot 4* 4th order laminar N-S 

PO IT-sat. urne-keU. 5c 
PO I T-sat u me- keO . 2 5c 
POlT-sat urne-ke0.25f 
POlT-saturne-rsinO.oc 
POIT-saturne-rsinO. 125c 

U. Poitiers 

( ,arpy, Manceau 

URANS, k-e model, 
dt.= 0.5, coarse grid 
URANS, k-e model, 
dt= 0.25, coarse grid 
URANS, k-e model, 
dt=0.25, fine grid 
URANS, RSM model, 
dt=0.5, coarse grid 
URANS. RSM model, 
dt— 0.125, coarse grid 

WARWK -K-neat-ke 
\VA RW IKK neat,- kenon 
WARWK ’K-neat-easm 

U. Warwick V U. Wales 

Preece, Tucker, Liu 

URANS, k-e model 
URANS, nonlinear k-e model 
URANS, EASM model 

VVASHU-wind-sa 

WASHU-wind-sst 

WASHU-wiud-sstles 

Washington U. 

Cui, Agarwal 

URANS, SA model 
URANS, SST model 
SST-LES 

NASA-t,lns3d-sa 

NASA-tlns3d-sa(coar) 

NAS A- 1, 1 ns3d-sa( fine ) 

N A S A- 1 1 ns3d-sa( low-d t, ) 

NASA-t.lns3d-sst 

NASA LaRC 

Vatsa, Turkel 

URANS, SA model 
URANS. SA model, 
coarse grid 
URANS, SA model, 
fine grid 

URANS, SA model, 
with lower dt 
URANS, SST model 
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Table 2. Case 1 Summary of Grids and Time Steps 


Label 

Grid type 

Grid size 

Time steps / cycle 

ONERA-flu3ni-3d-les 

3-D Structured, periodic 

930,000 cells 

5000 

O N E R A- fl u 3 m- 1 am 

2-D Structured 

51,700 cells 

5000 

0NERA-flu3m-sa 

2-D Structured 

51 ,700 cells 

5000 

UKY-ghost-sst 

2-D Structured Grid#l 

(53,553 points 

2880 

l ' I\ Y-ghost -sst ( f i lie ) 

2-D Structured Grid#2 

198,545 points 

2880 

GWU-vicar3d-3d(fine) 

3-D Structured, periodic 

(596,960 points 

14,000 

G W1 Lyjpai*3d-3d 

3-D Structured, periodic 

aha a a n — 

14,000 

GWU-vicar3d-3d(4.5d) 

3-D Structured, periodic 

696,960 points 

14,000 

NCAT-quasld+rans 

2-D Structured 

98,379 points 

118,567 

POIT-saturne-ke0.5c 

2-D Structured (no cavity) 

15,707 cells 

720 

P01T-saturne-ke0.25c 

2-D Structured (no cavity) 

15,707 cells 

1440 

POIT-saturne-keO .25f 

2-D Structured (no cavity) 

62,828 cells 

1440 

POIT-satu me- rsm 0 . 5c 

2-D Structured (no cavity) 

15.707 cells 

720 

POIT-saturne-rsm0.125c 

2-D Structured (no cavity) 

15.707 cells 

2880 

WA RWlCK-neat-ke 

2-D Structured (no cavity) 

4851 points 

3600 

WARWICK-neat-kenon 

2-D Structured (no cavity) 

4851 points 

3600 

WA RW I CK-neat-easm 

2-D Structured (no cavity) 

4851 points 

3600 

WASHU-wind-sa 

2-D Structured 

35,986 points 

10,000 

WASHU- wind-sst 

2-D Structured 

35,986 points 

10,000 

WA SHU- w i n d-ss t 1 es 

2-D Structured 

35,986 points 

10,000 

NASA-tIns3d-sa 

2-D Structured Grid#l 

63,553 points 

72 

NASA-t.lns3d-sa(coar) 

2-D Structured 

16,107 points 

72 

NASA-tliis3d-sa(fine) 

2-D Structured 

87,753 points 

72 

N ASA-tlns3d-sa(low-dt ) 

2-D Structured Grid#l 

63,553 points 

144 

NASA-tlns3d-sst 

2-D Structured Grid#l 

63,553 points 

72 
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Table 3. Case 2 Summary of Submissions 


Label 

Organization 

Authors 

Method 

NASA-cfl3d-sa 
N AS A-cfl3d-sa ( fine ) 
NASA-cfl3d-sst 
NASA-(’fl3d-easmko 

NASA LaRC 

Rumsey 

URANS, SA model 
URANS, SA model, fine grid 
URANS, SST model 
URANS, EASM-ko model 

llSTO-rans-t.lv 

USTO-raus-easm 

USTO fc ETH 

Azzi, Lakelial 

URANS, 2-layer k-e TLV model 
URANS, EASM model 

ON ERA-fl u3m-les 

ON ERA 

Dandois, Gamier, Sa.gaut 

LES 

( .'IRA-zen-ke-nocav 

C-IRA 

Marongiu 

URANS, k-e model 

NASA-fuii3d-sa 
N A S A- fu n3d-sa( fine ) 

NASA LaRC 

Atkins 

URANS, SA model 
URANS, SA model, fine grid 


Table 4. Case 2 Summary of Grids and Time Steps 


Label 

Grid type 

Grid size 

Time steps / cycle 

NASA-cfl3d-sa 

Full-plane, every other 
point of Structured Grid#l 

0.49 million cells 

720 

NASA-cfl3d-sa(fine) 

Full-plane, Structured 
Grid#l 

3.9 million cells 

720 

NASA-cfl3d-sst 

Full-plane, every other 
point of Structured Grid#l 

0.49 million cells 

1440 

NASA-cfl3d-easmko 

Full- plane, every other 
point of Structured Grid#l 

0.49 million cells 

1440 

USTO-rans-tlv 

Half plane structured 

0.21 million cells 

360 

L T STO-r ans-e as m 

Half plane structured 

0.21 million cells 

360 

ONERA-flu3m-les 

Full plane structured 

1.7 million cells 

6667 

CIRA-zen-ke-noeav 

Half plane structured 
(no cavity) 

775,680 cells 

720 

NASA-fun3d-sa 

Half plane unstructured 

46,000 nodes 

720 

NASA-fun3d-sa(fine) 

Half plane Unstructured 
Grid#! 

0.26 million nodes 

720 


14 of 31 


American Institute of Aeronautics and Astronautics Paper 2004-2217 


Table 5. Case 3 Summary of Submissions 


Label Organization 

A F RL-fd 1 3di-ke-c2 OAI fc AFRL 

A F R L-fdI3di-ke-c2-3d 


A FR L-fdl3di-ke-c4 


A FR L-fdl3di-ke-c4-3d 
A F R L-fd 1 3d i-ke-f2 


AZ-cobalt-des- l-3d 
AZ-cobalt-sa-l-bler 


Arizona State k 
Cobalt Solutions 


Authors 

Morgan, Rizzetta, 
Visbal 


Krishnan, Squires, 
Forsythe 


A Z-cobal t-sa- 1 - bl r 

A Z-coba l t-sa-l-r 

AZ-cobalt-sa-blc-3d 

AZ-eobalt-sa-blf-3d 

AZ-cobalt -sa-blni-3d 

AZ-cobalt-sst-blm-3d 

AZ-cobalt-sa-1 
AZ-cobalt -sa-lc 

AZ-cobalt -sa-2 
A Z-cobal t-sst-1 
A Z-cobal t-sst-lc 

AZ-cobalt-sst-2 
AZ-cobalt -sa'-ul 


A Z-cobal t-sa-u2 

BOEING-overflow-sa Boeing 

BOEING-overflow-sst 
M ETA-cfd ++ ke3-3d Metacomp Tech 
M ETA-cfd-h- hlns-3d 


Shmilovich, Yadlin, 
Clark 

SharifT, Batten 


Method 

BANS, k-e model, ‘2nd order 
RANS, k-e model, 2nd order, 3-D 
R.ANS, k-e model, 4th order 
RANS, k-e model, 4th order, 3-D 
RANS, k-e model, fine 
grid, 2nd order 

— dbJd 


RANS, SA model, B.L. on 
top wall, extended grid \v 
refined normal 
RANS. SA model, B.L. on 
top wall, grid vv refined normal 
RANS. SA model, grid w refined 
normal 

RANS, SA model, B.L. on 
all walls, coarse grid, 3-D 
RANS, SA model, B.L. on 
all walls, fine grid, 3-D 
RANS, SA model, B.L. on 
all walls, medium grid, 3-D 
RANS, SST model, B.L. on 
all walls, medium grid, 3-D 
RANS, SA model, fine grid 
RANS. SA model, coarsened 
normal direction 
RANS. SA model 
RANS, SST model, fine grid 
RANS. SST model, coarsened 
normal direction 
RANS, SST model 
RANS, SA model, unstructured 
fine grid 

RANS, SA model, unstructured 
~ ~ RANS, SA model 
RANS. SST model 
RANS, cubic k-e model, 3-D 
LNS. 3-D 
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Table 5. Case 3 Summary of Submissions, eont*d 


Label 

Organization 

Authors 

Method 

CTR-fluent-sa-1 

cm A fc CTR 

M arongiu , laccarino, 

RANS, SA model, fine grid 

( TR-fluent-sa-2 
( '1 R A-zen-sst-kprof-3 

( TRA-zen-sst-nj-3 

CIRA-zen-sst-oj-3 

CIR A-zen-sst-3 
CIR A-zen-ke-3 
CIRA-zen-kc 
CIRA-zen-sst 


( ’a tala no, Amato 

RANS, SA model 
RANS, SST model with 
k-profile at inflow 
RANS, SST model, 
normal-jet suction 
RANS, SST model, 
oblique- jet. suction 
RANS, SST model 
RANS, k-e model 
RANS, k-e model, fine grid 
RANS, SST model, fine grid 

N ASA-efl 3 cl- e as li i - 2 
NASA-cfl3d-sa-l 
NASA-cfl3d-sa-2 
N A S A - cfl 3 d-sst-2 
N A S A - cfl 3d - s a- t o p- 5 

NASA LaRC 

Runisey 

RANS. EASM model 
RANS. SA model, fine grid 
RANS, SA model 
RANS, SST model 
RANS, SA model, top wall 
modified for blockage 

N ASA-rans-sst-weno 5 

NASA LaRC 

Bala ku mar 

RANS, SST model, 
5th order weno 

NASA-fun2d-sa-ul 
N A S A - f u n 2d- sa - u 2 

NASA LaRC 

Vi ken 

RANS, SA model, fine grid 
RANS, SA model 

UAZ-cfl3d-easmfsm-3d 

U. Arizona 

Israel, Fasel 

FSM with linear EASM, 3-D 

UAZ-dns-3d 

U. Arizona 

Postl, Wernz, Fasel 

DNS, 3-D 

UK-ghost-sst-1 

UK-ghost-sst-2 

l . Kentucky 

Katani, Chen, Huang, 
Parimi. LeBeau, Huang 

RANS, SST model, 
fine grid 

RANS, SST model 

U M D-rans-sa-ceiit-2 
UMD-rans-sa-cent-2-3d 
UMD-rans-sst-cent-2 
U M D-rans-sa-cent- 1 

1 1 M D-rans-sa-roe-2 
I fMD-rans-sst-roe-2 

U. Maryland 

Duraisamy, Baeder 

RANS (central), SA model 
RANS (central), SA model, 3-D 
RANS (central), SST model 
RANS (central), SA model, 
fine grid 

RANS (Roe), SA model 
RANS (Roe), SST model 

US-fluent.-ko 
US-fluent-ke 
US-fluent-sa 
US-fluent-sst 
US-fluent.- v2f 

Utah State 

Spall, Phillips, Alley 

RANS, k-o model 
RANS, k-e model 
RANS, SA model 
RANS, SST model 
R ANS, v2f model 
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Table 6. Case 3 Summary of Grids and Time Steps 


Label 

Grid type 

Grid size 

Time steps / cycle 

AFRL-fdl5di-ke-c2 

2-D Structured 

50,410 points 

5,948 

A F R L- fd 13d i- ke-c2- 2d 

3-D Structured 
(half span / in viscid side) 

2.6 million points 

N/A 

A F R L- fdl :kl i- ke-c4 

2-D Structured 

50,410 points 

11,896 

AFRL-fdl3di-ke-c4-3d 

3-D Structured 
(half span / inviscid side) 

2.6 million points 

N/A 

AFR L-fdlMi-ke-f2 

2-D Structured 

199,790 points 

11,896 

A Z-cobal t-des- l-3d 

3-D Structured 
(periodic span) 

4.7 million nodes / 
4.5 million cells 

N/A 

AZ-cobalt-sa- 1-bler 

2-D Structured (extended 
upstream, refined 
normal dir) 

254.208 cells 

N/A 

AZ-cobalt-sa- 1-blr 

2-D Structured (refined 
normal dir) 

241,920 cells 

N/A 

AZ-cobalt-sa- 1-r 

2-D Structured (refined 
normal dir) 

241,920 cells 

N/A 

A Z-coba lt-sa-blc-3d 

3-D Unstructured 
(half span w plate) 

610,881 nodes / 
2.6 million cells 

N/A 

AZ-cobalt-sa-blf-3d 

3-D Unstructured 
(half span w plate) 

2.3 million nodes / 
10.7 million ceils 

N/A 

A Z-coba 1 t-sa- blm- 3d 

3-D Unstructured 
(half span w plate) 

1.1 million nodes / 
4.9 million cells 

N/A 

AZ-cobalt-sst-blm-3d 

3-D Unstructured 
(half span w plate) 

1.1 million nodes / 
4.9 million cells 

N/A 

AZ-cobalt-sa- 1 

2-D Structured Grid#l 

208.320 cells 

N/A 

AZ-cobalt-sa- lc 

2-D Structured 

1 10,880 cells 

722 

AZ-cobalt-sa-2 

2-D Structured Grid #2 

52,080 cells 

N/A 

AZ-cobalt-sst- 1 

2-D Structured Grid#l 

208,320 cells 

N/A 

AZ-cobalt-sst-lc 

2-D Structured 

110,880 cells 

N/A 

AZ-cobalt-sst-2 

2-D Structured Grid #2 

52,080 cells 

N/A 

AZ-cobalt-sa- ul 

2-D Unstructured 
Grid#l 

123,703 nodes / 
247,404 cells 

N/A 

A Z-coba 1 t-sa- u 2 

2-D Unstructured 
Grid#2 

57.152 nodes / 
114.302 cells 

N/A 

BOEING-overflow-sa 

2-D Structured overset 

47,790 points 

800 

B OEIN G-overflow-sst, 

2-D Structured overset 

47,790 points 

N/A 

M ETA-cfd + +ke3-3d 

3-D Structured 
(half span) 

2,472,520 cells 

N/A 

M ETA-cfd + +lns-3d 

3-D Structured 
(periodic span) 

1,816,000 cells 

N/A 
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Table 6. Case 3 Summary of Grids and Time Steps (cont'd) 


Label 

Grid type 

Grid size 

Time steps / cycle 

CTR-fluent-sa-l 

2-D Structured Grid#! 
w/o zone 4 

197,952 cells 

N/A 

(TR-fluent-sa-2 

2-D Structured Grid#2 
w/o zone 4 

49,488 cells 

N/A 

( -I RA-zen-sst-kprof-3 

2-D Structured Grid #3 
(no cavity) 

143,6)3 points 

N/A 

( -IR.A-zen-sst-nj-3 

2-D Structured Grid#3 
(no cavity) 

143,613 points 

N/A 

CIRA-zen-sst-oj-3 

2-D Structured Grid #3 
(no cavity) 

143,613 points 

N/A 

( 'IR A-zen-sst-3 

2-D Structured Grid #3 
(no cavity) 

143,613 points 

N/A 

CIRA-zen-ke-3 

2-D Structured Grid #3 
(no cavity) 

143,613 points 

720 

CIRA-zen-ke 

2-D Structured grid 
(no cavity) 

171,072 cells 

N/A 

CIRA-zen-sst 

2-D Structured grid 
(no cavity) 

171,072 cells 

N/A 

NASA-cfl3d-easm-2 

2-D Structured Grid#2 

52,080 cells 

N/A 

NASA-cfl3d-sa-l 

2-D Structured Grid#l 

208,320 cells 

360 

N AS A-cfl3d-sa-2 

2-D Structured G rid #2 

52,080 cells 

N/A 

N AS A-cfl3d-sst-2 

2-D Structured Grid#2 

52,080 cells 

N/A 

N A S A-cfl3d-sa- top-5 

2-D Structured Grid #5 

52.080 cells 

N/A 

N ASA- ran s- ss t- we n o5 

2-D Structured overset, 
(no cavity) 

103,452 points 

N/A 

NASA-fun2d-sa-ul 

2-D Unstructured 
Grid#l 

123,703 nodes / 
247.404 cells 

200 k 400 

NASA-fun2d-sa.-u2 

2-D Unstructured 
Grid# 2 

57,152 nodes / 
114,302 cells 

N/A 

UAZ-cfl3d-easmfsm-3d 

3-D Structured 
(periodic span, no cavity) 

2.8 million points 

200 

UAZ-dns-3d 

3-D Structured 
(periodic span, no cavity) 

105.2 million points 

N/A 

UK-ghost-sst-1 

2-D Structured Grid#l 
w/o zone 4 

197,952 cells 

N/A 

UK-ghost-sst-2 

2-D Structured Grid#2 
w/o zone 4 

49,488 cells 

N/A 
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Table 6. Case 3 Summary of Grids and Time Steps (cont'd) 


Label 

(j rid type 

Grid size 

Time steps / cycle 

UM D-rans-sa-cent-2 
I T M D-rans-sa-cent-2-3d 

U M D-rans-sst-eent-2 
V M D-rans-sa-cent- 1 
l 1 M D- rails- s a- roe- 2 
UMD-rans-sst- roe-2 

2-D Structured Grid#2 
3-D Structured Grid 
(half span) 

2-D Structured Grid #2 
2-D Structured Grid#l 
2-D Structured Grid#2 
2-D Structured Grid#2 

52.080 cells 
2.1 million cells 

52.080 cells 
208.320 cells 
52,080 cells 
52,080 cells 

012 

N/A 

N/A 

N/A 

012 

N/A 

US-fluent-ko 

2-D Unstructured Grid 

85.760 cells 

N/A 

US-fluent.-ke 

2-D Unstructured Grid 

85.760 cells 

N/A 

l T S- fluent -sa. 

2-D Unstructured Grid 

85,760 cells 

N/A 

US-fluent-sst 

2-D Unstructured Grid 

85.760 cells 

N/A 

US-fluent-v2f 

2-D Unstructured Grid 

85,760 cells 

N/A 
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Figure 1. Sketch of the case 1 comparison locations. 


Example determination of )et width 



Figure 2. Sketch showing definition of jet width for 

case 1. 
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Figure 3. Case 1 time histories of v- velocity at x=0. 
y=0.1mm. 
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Figure 4. Case 1 average v- velocity profiles at x=0. 




Figure 5. Case 1 average v- velocity profiles at 
y— 0.1mm. 
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v, m/s v, m/s 






Figure 6. 
y=4mni. 


Case 1 average v- velocity profiles at Figure 7 c ase i j e t width based on average v- velocity. 
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v, m/s V, m/s 


PIV Contours of V (135 deg) 



x (mm) 


WASHU-wtnd-sa 

WASHU-wind-sst 



Figure 10. Turbulence model effect on case 1 average 
v-velocity profiles at x=0. 

Figure 11. Case 1 contours of v-velocity, phaso=13 
(a) PIV experiment, (b), NASA-tlns3d-sa(fine). 
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Figure 14. Case 2 average u-velocity at y=0; (a) ID Figure 15. Case 2 average w-velocity; (a) z — lOmni ID 
downstream, (b) 2D downstream, (c) 8D downstream, downstream, (b) z=0.4inm and y=0, (c) z=10mm and 

y=o. 
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30 


phase=000 



Figure 16. Grid effect on case 2 u- velocity time histo- 
ries at x=50.63nim. y=0. z=0.4mm. 



Figure 17. Turbulence model effect on case 2 u- 
velocity time histories at x=63.5mm. y=0. z=10mm. 
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Figure 18. Case 2 phase- averaged u- velocity at y=0 
and ID downstream: (a) phase=0°, (b) phase=120°. 
(c) phase=240°. 
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phase=000 

NASA-cfl3d-sa 

NASA-cfl3d-sst 

NASA-cf13d-easmko 

USTO-rans-ttv 
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Figure 19. Turbulence model effect oil case 2 phase- 
averaged v* w* turbulent stress at x=44.45mm. y=0. 
phase=0°. 



Figure 20. Example contours of phase- averaged u- 
velocity. ID downstream. phase=120°; (a) PIV exper- 
iment, (b) NASA-cfl3d-sa(fine). 
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Figure 21- Sketch of the case 3 comparison locations. 
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Figure 22. Summary of case 3 results compared to Figure 23. Summary of case 3 results compared to 
experimental data for no-flow-control condition; (a) experimental data for suction condition: (a) C'p, (b) 
C'p, (b) separation location, (c) reattachment location. separation location, (c) reattachment location. 
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Figure 24. Surface pressure coefficients for case 3 re- 
sults that accounted for tunnel blockage. 


Suction, x/c=1.0 



No- flow- control. x/c=1 .2 



Suction, x/c=1 .0 



Figure 26. Sample of case 3 suction condition u- 
velocity profiles at x/c=1.0 (downstream of experi- 
mental reattachment ) . 


Figure 25. Sample of case 3 no-flow-control condition 
u-velocity profiles at x/c=1.2 (downstream of experi- 
mental reattachment). 
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Suction, xJc=Q& 



u/Uinf 


Suction, x/c=0.8 




Figure 28. Summary of case 3 long-time-average Cp 
results for oscillatory control condition. 


Figure 27. Sample of case 3 suction condition results 
at x/c=0.8 (inside separation bubble); (a) u-velocity, 
(b) turbulent shear stress. 
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